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The chirally-odd twist-3 parton distribution function e'^{x) of the nucleon is studied in the large- 
Nc limit in the framework of the chiral quark-soliton model. It is demonstrated that in spite of 
properties not shared by other distribution functions, namely the appearance of a 5(a;)-singularity 
and quadratic divergences in e'^(x), an equally reliable calculation is possible. Among the most 
remarkable results obtained in this work is the fact that the coefficient of the (5(a;)-singularity can be 
computed exactly in this model, avoiding involved numerics. Our results complete existing studies 
in literature. 

PACS numbers: 12.39.Fe, ll.15.Pg, 12.38.Lg, 13.60.Hb 

Keywords: chirally odd twist-3 parton distribution, large Nc limit, chiral soliton model of the nucleon 

I. INTRODUCTION 

The twist-3 parton distribution function e°(x) [1] was so far subject to modest interest in hterature, in spite of its 
remarkable theoretical properties, because of its chirally odd nature which makes it difficult to access in experiment. 
The probably most striking theoretical property is the existence of a i5-function-type singularity at a; = in e°'{x) 
which follows from the QCD equations of motion [2-4]. 

Since it is chirally odd e°'{x) can contribute to an observable only in connection with another chirally odd object. 
For example, e"(x) and the chirally odd Collins fragmentation function [5] contribute to an azimuthal asymmetry 
in semi-inclusive deeply inelastic scattering of longitudinally polarized electrons off unpolarized nucleons [6, 7] which 
was measured [8-10] and used to extract first information on e°'{x) [11]. Later it became clear that e"(x) and H^^ 
are not the only contribution to this observable [12-16]. Studies of two-hadron production in semi- inclusive deeply 
inelastic scattering could provide a more direct and easier access to e°'{x) [17]. 

In experiment the (5(x)-singularity could be observed only indirectly, namely as a discrepancy for the first moment 
of e°'{x) between the theoretical result, which includes the point x = 0, and an experimental result where only x > 
contribute [4]. A "direct observation" of the (5(a;)-contribution, however, is possible in models. 

The appearance of a 5(a:)-singularity makes e°'(x) a particularly interesting but, at the same time, also demanding 
object for model studies. 

While there is no such singularity in bag [1, 18] or spectator [19] models, a (5(a;)-contribution in e°'{x) was found 
in (1 -I- l)-dimensional toy model or perturbative one-loop calculations [20, 21], and in the chiral quark-soliton model 
(xQSM) [22-24]. The latter is the ground for the present study. 

The existence of a (5(a;)-contribution in e°'{x) in the xQSM was proven independently in [22, 23]. Moreover in [22] 
a first approximate calculation of e°(x) was presented and a method was suggested, how to compute in practice in 
this model a parton distribution function in which a 5(a;)-singularity appears. This method was explored to calculate 
e''(x) exactly in the xQSM in [24]. 

The calculation of e°'{x) in the xQSM is demanding not only because of the appearance of the (5(a;)-contribution. A 
further complication is due to the fact that e"(x) is quadratically UV-divergent, in contrast to the parton distribution 
functions studied so far which are UV-finite or at worst logarithmically divergent. For these reasons it is worthwhile 
to present an independent exact calculation of e°'{x) in the xQSM. This is the purpose of the present study. 

In order to supplement and complete previous works [22-24] we shall devote particular effort to the demonstration 
that the involved numerics in this calculation is well under control. The tool we shall use for that is the so-called 
gradient expansion. This work, however, goes beyond a mere reexamination and confirmation of the previous studies 
[22, 24]. In fact, among the new insights, the most remarkable one is the observation that the result for the coefficient 
of the S{x)-term in e°(x), obtained by means of the gradient expansion in [22], is exact. 

This note is organized as follows. In Sec. II we review the general properties of e°'{x). In Sec. Ill we introduce the 
xQSM, and discuss how e°'{x) is described in this model. The Sees. IV, V and VI are devoted to the singlet fiavour 
combination which is leading in the large- iVc hmit, and we discuss in detail respectively its singular and regular parts, 
and the origin of the 5(a;)-singularity. In Sec. VII we study the non-singlet flavour combination which is the subleading 
structure in the large- A'c expansion. In Sec. VIII we summarize and discuss our results which we compare to previous 
studies in the xQSM and other models of e°-{x) in Sees. IX and X, respectively. Sec. XI contains the conclusions. 
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II. THE DISTRIBUTION FUNCTION e"(a;) 

Let us first recall some general theoretical properties of e°'{x), a more extensive review can be found in Ref. [4]. 
The twist-3 distribution functions e'(a;) for quarks of flavour q and e^{x) for antiquarks of flavour q are defined as [1] 

with the light-like vector n and [0, An] denoting the gauge-link. The dependence on the renormalization scale fi, which 
we do not indicate in (1) for brevity, was studied in [2, 3]. 

The equations of motion of QCD allow to decompose e*(a;) as [2-4] (in this work we shall consider e'(a;) in the 
chiral limit, and limit ourselves to only indicate current quark mass effects) 



e'(^) = Wtr {N\MO)Mm) + eLsix) +0{^]. (2) 



S{x) 

2Mn ' -twj^-/ ■ - yj^^ 

The contribution ej^^g(a;) in Eq. (2) is a quark-gluon-quark correlation function, i.e. the actual "pure" twist-3 ("in- 
teraction dependent") contribution to e°'{x), and has a partonic interpretation as an interference between scattering 
from a coherent quark-gluon pair and from a single quark [1]. Its first two moments vanish. Consequently, the first 
moment of e'^{x) is saturated by the ^(a;)-contribution and the second is due to quark mass effects only 



\ '^'^ ^'^''^ = 2]^^^l ^5,(0)^,(0) \N) , (3) 

777 

dxxeoix) = j^N,. (4) 



In Eq. (4) Nq denotes the number of the respective valence quarks (for proton iV„ = 2 and N^, = 1). 

The first moment (3) of the flavour singlet combination is related to the pion-nucleon sigma-term ct^n as 



dx (e« + e'')(x) = ^ = ^ 
m m 



(5 - 15) , (5) 



where m = ^ {niu + m^) and a{t) denotes the scalar isoscalar nucleon form-factor which is known only at the Cheng- 
Dashen point t = 2rn^ where it is related by low-energy theorems to pion-nucleon scattering amplitudes. Analyses 
yield (T(2m^) = (56 — 88) McV [25, 26]. The difference (7(2m^) — (t(0) = 14MeV was obtained from a dispersion relation 
analysis [27] and a similar result was found in chiral perturbation theory calculations [28]. Thus a^^ ~ (42 — 74) MeV 
and with m = (7 ± 2) MeV [29] one obtains the result quoted in Eq. (5) which refers to the scale = 1 GeV. 
The flavour non-singlet combination satisfies an equally interesting sum rule, namely 

' dx (e" - e'^) (X) = (^"-^P)hadr ^ _ 

.1 md - ruu 

Here (M„ — Mp)hadr ~ 2 MeV [29] denotes the hadronic mass difference between neutron and proton, i.e. the mass 
difference in the absence of electro-weak interactions. With — m„ « 5 MeV [29] one obtains the result quoted in 
Eq. (6) at a scale of /x = 1 GeV. It is important to keep in mind that the numbers in Eqs. (5, 6) are not due to the 

"valence" structure of e'^(.x), but solely due to the (5(.T)-contributions. 
In the large- iVc limit the behaviour of the different flavour combinations of e"(a;) is as follows [4, 22] 

(e" + e'')(x) = Nld+{N,x) , 

(e"-e^)(x) = N, d-{N,x) , (7) 

where d±{y) are stable functions in the limit Nc ^ oo for fixed arguments y = NcX. This implies the following 
hierarchy 

l(e" + e'^)(x)|»l(e"-e'^)(x)l (8) 

in the large- A'^c limit. Notice that the twist-2 distribution function /f (cc) exhibits the analog flavour dependence. 
Noteworthy, e'(a;) and ff{x) become equal in the non-relativistic limit [4, 22]. 
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III. CHIRAL QUARK SOLITON MODEL, ITS APPLICATIONS, LIMITATIONS, AND e"(a;) 

The effective theory underlying the chiral quark soUton model (xQSM) [30] was derived from the instanton model 
of the QCD vacuum [31, 32], and is given by the partition function [33] 

(9) 



■gff = VtpVip VU exp 



i <rxti){i^ -MV"- -m)'4) 



where ~ cxp(?'75T"7r"). M denotes the dynamical quark mass due to the spontaneous breakdown of chiral 
symmetry, and m is the current quark mass responsible for explicit chiral symmetry breaking effects. The effective 
theory (9) contains the Wess-Zumino term and the four- derivative Gasser-Leutwyler terms with correct coefficients. 
This effective theory provides an approximation to the dynamics of light quarks valid at low energies below 

p-^Ri600MeV, (10) 

where Pav is the average instanton size. Corrections to this picture are of the order (Mpav)^ ~ 30%. It is important 
to remark that (Mpav)^ is proportional to the parametrically small instanton packing fraction 

(Mpav)' (X « 1 , (11) 

where i?av denotes the average distance between instantons. The numerical smallness of the parameter Pav/Rav ^ 3 
played an important role in the derivation of the effective theory (9) from the instanton vacuum model [31]. 

Let us briefly recall how the effective theory (9) describes nucleons [30] . In the leading order of the large- A^c limit 
the pion field is static and the soliton energy E'soi is a functional of this field given by 



E„<0 



(12) 

reg 



The En in (12) are the eigen-energics of the onc-particlc Hamiltonian 

H\n) = En\n) , H = -i'y°j''dk + -(°MW + 7°m , (13) 

and the Ena of the free Hamiltonian which follows from (13) by replacing U'^-'' 1. The spectrum of the Hamiltonian 
(13) consists of an upper and a lower Dirac continuum, and — for a strong enough pion field — of a discrete bound 
state level of energy E\f,^. By occupying the discrete level and the states of the lower continuum each by quarks 
in an anti-symmetric coknir state and subtracting the vacuum, one obtains a state with unity baryon number. The 
quantity i?soi[C^] is logaritliinically divergent and has to be regularized as indicated in (12). The minimization of 
-£'soi[t^] is performed for symmetry reasons in the so-called hedgehog Ansatz 7r°(x) = e" P{r) and determines the 
self-consistent soliton profile Pc{r) and soliton field Uc- The nucleon mass Mn is given by Eso\[Uc]- 

Nucleon states with a definite momentum and quantum numbers are obtained by considering translational and 
rotational zero modes of the soliton. In order to include corrections in the 1/A^c-cxpansion one considers time- 
dependent pion field fiuctuations around the self-consistent solution. Hereby one restricts oneself to time-dependent 
rotations ?7c(x) R{t)Uc{x)R^t), where the collective coordinate R{t) is a rotation matrix in SU(2)-flavour space. 
In this approximation the integral over all pion field fluctuations in (9) is given by the path integral over the collective 
coordinates and solved to leading order in the collective angular velocity Q = —iR^dtR. The expansion in n is 
justified, since the corresponding soliton moment of inertia 

6 Ei-En ^ ' 

J, non ■> 

is of 0{Nc) and thus large, such that the soliton rotates slowly. The sums in (14) go over occupied "occ" states n (non- 
occupied "non" states j), i.e. over states with En < Ei^v {Ej > Eicv)- The xQSM provides a particular realization of 
the general large- A^c picture of the nucleon [34]. 

With M = 350 MeV fixed from instanton phenomenology [31] and the precise value of the cutoff (10) adjusted 
within the chosen regularization scheme to reproduce the experimental value of the pion decay constant, there are 
no free parameters in the xQSM. In this sense the model allows to evaluate in a parameter-free way nucleon matrix 
elements of QCD quark bilinear operators (A'^|'i/'(2^i)r^/'(2:2)lA'^) where F is some Dirac- and fiavour-matrix. In this way 
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numerous baryon properties such as electromagnetic, axial or scalar form- factors, etc., were computed in the model 
and found to agree with data to within an accuracy of typically (10-30)% [35 41]. (Notice that in some studies M was 
varied, and e.g. for M = 420 MeV in the proper-time regularization some static properties were found to be better 
described [35]. However, the dependence of the results on variations of M is moderate, and within model accuracy.) 
An interesting recent application includes the extension of the model to the description of the implicit dependence of 
the nucleon mass on the pion mass in the regime studied in lattice QCD [42, 43]. 

The model was also applied to studies of usual [44-50] and generalized [51-54] twist-2 quark- and antiquark- 
distribution functions at a low scale of 0{p~^). The small parameter in Eq. (11) is of crucial importance for justifying 
the computation of twist-2 (generalized) parton distribution functions in the xQSM [32]. To leading order of the 
instanton packing fraction (11) the model quarks can be identified with QCD quark degrees of freedom, while gluon 
degrees of freedom appear suppressed by this parameter [32]. This is essential to guarantee a consistent description, 
and the model results satisfy all general QCD requirements such as sum rules, inequalities, polynomiality [44-54] and 
agree — as far as those quantities are known — to within (10-30)% with parameterizations [55, 56]. 

Is the model also applicable to studies of twist-3 distribution functions? The answer to this question cannot be 
found within the model itself, since twist-3 quantities can be rewritten by means of QCD equations of motion in 
terms of contributions in which gluon fields appear explicitly. Instead it is necessary to consider this question in the 
instanton vacuum model, i.e. in the theory from the model is derived. 

The twist-3 distribution functions g^ix) and h1{x) were studied in the instanton vacuum model [57] and it was 
found that the pure twist-3 parts in the Wandzura-Wilczek(-Iike) decompositions of these functions are strongly 
suppressed by powers of the instanton packing fraction (11). As g^{x) and h1{x) can be defined without the explicit 
appearance of gluon fields — just as e°(x) in Eq. (1) — these functions can in principle be computed in models 
without gluon degrees of freedom [1]. In the xQSM this was done in [58] and the pure twist-3 parts of g^{x) and 
h1{x) were found to be small. Thus, in these cases the xQSM respects the results of the instanton vacuum model, 
i.e. of the theory from which it was derived. 

It would be interesting to know whether calculations of e"'{x) in the xQSM are similarly in agreement with the 
instanton vacuum model. Answering this questions is, however, beyond the scope of this work. Here we shall take the 
practical point of view of Ref. [1], and explore the special property of e"(x) that allows to define it in terms of quark 
fields only, without resorting explicitly to gluon degrees of freedom, which makes it possible to study this quantity in 
effective approaches with quark and antiquark degrees of freedom only. 

In the large-A^c limit different flavour combinations of nucleon quantities appear usually at different orders in 
the large- A''c counting. This is also the case in the xQSM, which respects all large- A^c counting rules. For e"(x) the 
isoscalar flavour combination is leading, and the isovector one appears only at subleading order in the 1 /Nc expansion. 
The model expressions for the different flavour combinations of e"'{x) in proton read [22] 

(e" + e'^)(x) = AfcMN^(nl7°5(:cMN-p^-^„)ln)reg (15) 

nocc 

= - A-cMn Y1 {nh°S{xM^ -f- £„)|n),eg (16) 

n non 



j all j'#n 

A^cMn ( 2 d 



127 . V -^7 - En dxMjs 



{n\T''\j){jKl"S{xM^-p^ - En)\n) (18) 



where vacuum subtractions analog to (12) are understood. 

The possibility of computing model expressions for parton distributions in the two different ways, by summing 
over occupied (15, 17) and non-occupied (16, 18) states is deeply connected to the analyticity properties of model 
expressions and founded on the locality of the model [44] . 

Therefore it is of importance to demonstrate explicitly that the equivalent formulae, (15, 16) and (17, 18), yield 
respectively the same results. This not only provides a very useful test of the numerics. The explicit demonstration 
of the "equivalence" of the different representations provides a crucial test for the internal, theoretical consistency of 
the model itself, and we shall devote much effort to this point. 



5 



IV. CALCULATION OF (e" + e'^){x) 



The isoscalar combination (e" + e'^){x) contains a (5(x)-singularity, as 
proven in [22] and independently shown in [23]. A practical procedure 
to cope numerically with such a singularity in the xQSM was suggested 
in [22], and used in [24] to confirm numerically the existence of the 6{x)- 
singularity. In this Section we present an independent study, which will 
confirm the findings of Ref. [24]. Hereby we shall focus on the demon- 
stration that the involved numerical calculation is well under control. 

The contribution of the discrete level to (e" + e'^){x) as well as to any 
quantity in the model UV-finite. It can be computed by directly solving 
the eigenvalue problem (13) for the discrete level state [44], or by using 
the method described below. The result is shown in Fig. 1. 

The continuum contribution can be computed in two ways, which fol- 
low from Eqs. (15, 16) (and to which we shall, for simplicity, continue 
referring as sums over occupied and non-occupied states), namely 



(e"+e'^)(x)i^^ 




(e"+e'*)(x)eo„t,reg = iV,MN ^ {72\j'> SixM^ - - EM^cg , (19) 

= -N,M^ (n|7°'^M^N-p3-^„)|n)„g . (20) 



-1 -0.5 0.5 

FIG. 1: The contribution of the discrete 
level to (e" -I- e'^){x) as function of x. 



E„>o 

It is quadratically UV-divergent and has to be regularized, as indicated in (15, 16) and (19, 20). The Pauli-Villars 
subtraction method is the only regularization known in the xQSM which preserves all general properties of distribution 
functions (QCD sum rules, positivity, etc.) [44]. 

In order to regularize the continuum contribution and to ensure the equivalence of the summations over occupied 
and non-occupied states in Eqs. (19, 20), see [22] for a detailed discussion, two Pauli-Villars subtractions are necessary 



with 



(e'' + e'^)(a;)cont,rcg - (e" + e'')(x, M^nt - ai(e" + e''){x, Mi)eo„t - a2(e" + e''){x, M2)c 

M Ml - ]VP 



ai 



M Af| 



a2 



Ml M| - Ml ' ' M2 M'i - M'l 



M2 > Ml > M 



(21) 



(22) 



The values for the Ali are fixed by regularizing the model expressions for the pion decay constant and the vacuum 
quark condensate (ipip) which are given in the effective theory (9) by the (Euclidean) loop integrals 



J TT 



(27r)4 ipl + M^y 



{il}^)) = (vacKV'tiV'" + V'd'0d)|vac) = 



dV 8iVcM 



(2^)4 pI + M2 



(23) 



The model expression for the pion decay constant is regularized and its phenomenological value /tt = 93MeV is 
reproduced by a single Pauli-Villars subtraction with Mi = 556 MeV. Two subtractions analog to (21, 22) are needed 
to regularize the model expression for the quark vacuum condensate (23). In order to reproduce the phenomenological 
value (V^V') = - (280 ± 30)3 MeV^ quoted in [29] one should use M2 = (2.11q 7) GeV. For reasons which we will explain 
below a small value of M2 is preferable, and we choose M2 = 986MeV which yields {'tp'4') = — (220MeV)'^. This is 
sufficiently close to the phenomenological value — considering the typical accuracy of the model. 

In order to evaluate {e'" + e'^){x) we use the following procedure [45]. For the numerical calculation the Hamiltonian 
(13) is expressed and diagonalized in the basis of the eigenstates of the free Hamiltonian. The spectrum of (13) 
is discretized by placing the soliton in a finite but sufficiently large spherical box of the size i?box and imposing 
appropriate boundary conditions [59]. The spectrum is made finite by cutting off quark momenta above some large 
numerical cutoff Anum chosen much larger than any other (physical or numerical) scale involved in the problem. We 
use i?box = (8 - 12) fm and A„ui„ = (8 - 9) GeV. 

To compute the continuum contribution wc introduce an intermediate regularization for any of the contributions 
in Eq. (21) schematically as 



-e^){x,M,) = N,M^ 



J2{n\j°S{xM^ -f- E,,)\n) R{E„, A) 



(24) 



Mi 
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where it is understood that the expression in squared brackets is to be evaluated with the Hamiltonian (13) if Mi = M, 
or with Hamiltonians analog to (13) but with M replaced by Mi or M2, and where vacuum subtraction is implied. 
In (24) R{lu, A) is a smooth regulator function with i?(0, A) = 1 and R{uj, A) ^ for \uj\ 00. The intermediate 
cutoff A must satisfy Mi ^ A <C Anum- It is due to this condition that we prefer a low value of M2, see above. In 
practice we use A ~ (3 — 6) GeV. The regulator function can be chosen to be of e.g. Gaussian or Wood-Saxon type. 
The dependence on A and the choice of R is removed at the end of the calculation by an extrapolation procedure. 

It is convenient to turn (24) into a spherically symmetric form by replacing — > u • p where u is a unit vector. 
(We recall that the 3-direction was singled out arbitrarily by choosing the spatial component of the light-like vector 
n in (1) along that axis.) Averaging over the possible orientations of u yields 



(e" + e'')(a;,M,) = iVcMN 



^©(IpI ' l^^-^N - En\)\n) R{En, A) 



(25) 



Mi 



As we work in a discrete basis the expression (25) is a discontinuous function of x due to the appearance of the 
G-function, and would become continuous only in the infinite volume limit. Rather than trying to take this limit 
numerically, which would be a time-consuming procedure, one may instead smear the expression in (25) by convoluting 
it with a narrow Gaussian as 



(e" + e'^)(a;,MO 



77r 



1/2 



dx' e-(^-^')'''" 



\x'MN~En\)\n) R{En, A) 



(26) 



Mi 



In the limit 7^0 one recovers the original expression (25). The parameter 7 > has to be chosen such that it is, on 
the one hand, sufficiently large compared to the typical splitting of energy levels in the discretized spectrum, and on 
the other hand, sufficiently small such that the "smeared function" is still a good approximation to the true result. 
For our box sizes 7 = 0.1 is adequate [45]. In the end of the day the smearing can be removed by a deconvolution 
procedure, though in practice one finds that continuous functions are sufficiently well approximated by (26). 

It is precisely the smearing procedure which enables one to cope numerically with a (5-function-type singularity. In 
fact, the smearing trick, turns the (5(a;)-contribution in (e^ + e'^)(a;) into a narrow-Gausssian of a well-defined width 
7 centered around x = 0. Using the self-consistent profile from [46] which yields 1140 MeV for the soliton energy 
and the above described parameters we obtain for the regularized continuum contribution (e" -I- e'^)(a;, A)cont for the 
intermediate cutoffs A = 3 and 6 GeV the results shown in Figs. 2a and b. 

The results depend on whether one sums over occupied, Eq. (19), or non-occupied, Eq. (20), states and on the 
intermediate cutoff A. However, after extrapolating A ^ 00 both formulae (19) and (20) yield within a numerical 
accuracy of 1% the same result, which is shown in Fig. 2c as solid line. 




FIG. 2: The continuum contribution to (e" + e'*)(x) as function of x. (a) The result from the sum over occupied levels (19) for 
different values of the intermediate cutoff A. (c) The same as (b) but from the sum over non-occupied levels (20). (d) Solid 
line: the final result obtained after the extrapolation A — > 00. Dotted line: the result from the gradient expansion [22], see text. 
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V. THE ORIGIN OF THE <J(a;)-SINGULARITY IN (e" + e'^){x) 

Having convinced ourselves that one obtains for the continuum contribution of (e" + e'^){x) the same result, irre- 
spective whether one computes it by means of (19) or (20), we have confidence in our numerical results, and are in 
the position to address the question: What precisely gives rise to the (5(a;)-singularity in (e" + e'^){x) in the xQSM? 

As there is apparently no S{x)-teTm in the level part, see Fig. 1, one has to focus on the continuum contribution. 
The continuum contribution (19) or (20), is given in the gradient expansion by [22] 

(e" + e'^) {xUnt = C 5{x) + 0{VU) , C = {iPi>) J (^cos P{r) - 1^ , (27) 

where 0{W) denotes terms which (i) contain one or more gradients of the chiral field U, i.e. are suppressed in a 
chiral expansion, and (ii) are regular functions of x. 

The result (27) is remarkable for two reasons. First, the coefficient C of the 5(a;)-function can be computed exactly. 
Second, the quadratic and logarithmic divergences appearing in C can be reexpressed in terms of (ipip) [22]. 

It is is interesting to confront (27) which describes the singular part of (e" + e"^) (x) in the gradient expansion exactly 
with our numerical result. For that we evaluate C using the same Pauli-Villars masses, i.e. {"tptp) = —(220 Me V)^, and 
the same self-consistent profile P{r) as in the numerical calculation which yields 

Chare = 4.88 . (28) 

(The introduced index "here" reminds that this result holds for the Pauli-Villars masses used here, in this calculation.) 
Moreover, we smear the (5(x)-function according to (26) with the same parameter 7 = 0.1 as used in the numerics. In 
this way we obtain the result shown as dashed line in Fig. 2c. 

The impressive agreement in Fig. 2c fully confirms the findings of Ref. [22] that the (5(.T)-contribution in the xQSM 
originates solely from the leading order of the gradient expansion where it is related to the quark vacuum condensate. 
Fig. 2c also shows that the regular contribution to the continuum part of (e"-|-e'^)(x) is small (on the scale in Fig. 2c). 
We observed similar agreements by varying the numerical parametcirs (Pauli-Villars masses Mi, profile P{r), 7, etc.). 
These observations illustrate the utility of the gradient expansion as a powerful tool to control the numerics. 

Let us report the following detail which further increases our faith into the quality of the numerics. In order to 
achieve the satisfactory agreement in Fig. 2c we considered the finite size of the spherical box used in the numerics, 
and integrated in (27) over the radial component r = |r| only up to r = i?box — 12 fm which gave the result in (28). 
In the chiral limit the profile function behaves as P{r) = A/r^ at large r, where A is related to the isovector axial 
coupling constant by ^ = {2>g a) / {^t^ f^) ■ In practice the asymptotics sets in already for r > 3fm [46]. Thus, the 
finite size effect for the coefficient C in (27) is 

SC = {ipi)) 47r / dr cos P(r) - 1 = - (V'V') ^ • (29) 

The finite size effect slowly vanishes with increasing box size. For -Rbox = 12 fm one has (5Chcro = 0.196, i.e. the result 
for Chere in (28) is about 4% smaller compared to its infinite volume limit. Neglecting this effect would yield a clearly 
visible mismatch in Fig. 2c. Thus, as a byproduct we see that finite box size effects in our calculation are of the order 
of magnitude of a few percent, which is acceptable considering the typical model accuracy. 

Thus, besides confirming numerically the presence of a (5(a;)-contribution in (e" + e'')(a;) as was done previously in 
[24], we furthermore have numerically confirmed the fact that the coefficient C of the 5(a;)-function can be computed 
exactly using gradient expansion [22] . This is of importance because only the gradient expansion allows to relate C 
to the quark vacuum condensate. 

The exact prediction for the coefficient C in the xQSM is therefore 

Cexact = lO-S^si (30) 

using the self-consistent profile [46] and the phenomenological value ('^V) = — (280±30)^ MeV^ [29] whose uncertainty 
yields to the error shown in (30). 

We remark that in [24], where the Pauli-Villars masses were fixed to reproduce {'fp'tj}) = — (280MeV)^, the result 
C ~ 9.9 was obtained numerically. This is in good agreement with the central value of (30) and the small discrepancy, 
apart from numerical uncertainties (finite box size effects), is due to the slightly different value of the constituent 
mass M = 375 MeV used in [24], which yields a somehow different self-consistent profile. We checked that using larger 
values of M2 ^1.5 GcV we are able to reproduce larger values of C. The price to pay for that is, however, a worse 
numerical stability because then the required hierarchy Mi ^ A Anum holds less satisfactorily, unless one increases 
A and Anum accordingly from which we refrained being limited by the computing resources available to us. 

Having established the presence of a (5(.T)-type-singularity in [e^ + e'^){x), confirming thereby the findings of [24], 
we now turn to the task of computing the regular part of the continuum contribution. 
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VI. CALCULATION OF a;(e" + e'^){x) 



Being interested in the regular part of the continuum contribution to (e" + e'^){x) the method of Sec. IV is not 
adequate. The smearing trick, which allows to visualize the delta-function, is of disadvantage in this case. It turns the 
(5(a;)-function into a Gaussian that penetrates into the regions x and dominates there over the regular part, Fig. 2c. 
Thus, a reliable computation of the regular part requires a different technic. Here we shall compute a;(e" + e'^)(a;) 
where the S{x) drops out. (Due to the smearing procedure (26) computing a;(e" + e'^)(a;) is, of course, not the same 
as computing (e" + e'^){x) and multiplying the result by x.) 

First we have to clarify how to regularize a;(e" + e'^){x). It is worth to reconsider this point because x{e^ + e'^){x) 
and (e'* + e'^)(a:) differ by the appearance of a 5{x)-terni. As shown in [22], the double subtraction (21, 22) is needed 
to make finite the coefRcient of the (5(a;)-function. But is it also necessary to regularize a;(e" + e'^)(x)? It is worth to 
reconsider this point: one could save a lot of computing time, if no or only one subtraction were necessary. 

By going step by step through the Eqs. (32-40) of Ref. [22], one believes, at a first impression, that a single 
subtraction (corresponding to (21, 22) for M2 — > 00) is enougii to remove from a;(e" + e'*)(a;) a quadratic divergence 
and restore the equivalence of summations over occupied and non-occupied levels [22]. However, such a conclusion is 
premature and we show here that a double subtraction is adequate. 

In order to illustrate why we integrate a;(e" + e'^)(a;) from (19, 20) over x, explore the model equations of motion, 
and arrive at expressions for the continuum contribution to the second moment of (e" + e'^) {x) which read [22] 

NcM 



I 



da;a;(e"-|-e'*)(a;)cont 



2Mi 



N 



J2{nm + U^)\n)r 



-En<0 

NcM 



2Mi 



N 



(n|(C/ + C/t)|n) 



reg 



E„>0 



Clearly, the difference of the two expressions in (31) and (32) must be zero, i.e. one expects 



ireg 



2Mn 



5](n|(f/ + f/t)|n) 



reg 



n all 



iVcM 

Mm 



Sp 



U + W 



(31) 



(32) 



(33) 



n) 



Here Sp denotes the functional trace which can be saturated by any complete set of functions, Sp[ . . . ] = X]n('^l 
being one example. Notice that the "—1" in Eq. (33) is due to the explicitly included vacuum subtraction. 

In the numerical calculation the expression for ylj-og is evaluated with an intermediate regularization, see Eq. (24), 
such that in the intermediate step we are interested in the following expression 



A(A.«) = ^Sp 



U + W 



R{H, A) - R{Ho, A) 



(34) 



We evaluate A{A,M) in Eq. (34) in gradient expansion where R{H, A) = R{Ho, A) -I- C(VC/) holds, and saturate 

the trace over Dirac- 
= i?(— w, A), we obtain: 

i?(\/k2 + M2,A) . (35) 



Sp[ - • •] = / (2^^3 (kjtr . . . jk) in the spectrum of the free momentum operator where tr denotes the trace over Dirac- 
and flavour-indices. As positive and negative energies must be regularized equally, i.e. R{co, A) 

8NrM f ,^ , _ /■ d-'^k 



A{A, M) 



d'^x (cos P(r) 



(27r) 



We recognize, most easily by employing a definite regulator such as e.g. a simple Gaussian i?(w. A) = exp(— w^/A^), 
that ^(A, M) has a cubic and a linear divergence for A ^ 00. The double subtraction (21, 22) not only removes these 
divergences but also restores the equivalence of the summations over occupied or non-occupied states, i.e. 



ireg 



lim A{K) = , where A(A) = A{A, M) - aiA{K, Mi) - a2^(A, M2) . 



(36) 



One can convince oneself similarly as done in [22] that higher orders in VC/ omitted in (35) do not spoil the above 
argumentation. Thus, a;(e"-|-e'')(a;) must be regularized exactly in the same way as (e" -|- e"^) (x) according to (21, 22). 

The reason why the calculation in Eqs. (32-40) of [22] lead us to the premature conclusion that a single subtraction 
could be sufficient is because here and in [22] different matrix elements were evaluated in gradient expansion, namely 



{n\En^''\n) 



used in [22] 



^2=- l(n|{ff,7°}|n) = ^M(n|(C/ + C/t)|n) 



used here 



(37) 



which are connected by equations of motion (eom). The latter are, of course, not respected in a truncated expansion. 
Another example that the gradient expansion yields results at variance with eom is documented in Sec. 7 of Ref. [44] . 
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Having clarified the issue of regularization we turn to the numerical 
evaluation. At first glance one may have the impression that different 
expressions exist, namely 



x(e"+e'')(x,A)cont 



n occ 



(38) 



E 

n occ 



{n\f6{xMn -p^- En) if + En) \n). 



and analog expressions with summations over non-occupied states. How- 
ever, it is gratifying to observe that after averaging over directions either 
of these expressions just yields x times the result in Eq. (25). 

The numerical procedure of Sec. IV yields for a;(e" -I- e'')(a;, A)cont the 
results shown in Fig. 3. We make two interesting observations. First, 
compared to the huge effect of the smeared out (5(x)-function in Fig. 2 
the regular contribution is small. Second, the results for the sums over 
positive and negative energy states exhibit a tendency with increasing 
intermediate cutoff A to converge slowly towards a common result. How- 
ever, this convergence is remarkably slow, and the numerical stability of 
a A ^ oo extrapolation is poor, especially for small < 0.2. 

The slow convergence can be understood as follows. We compute 
the moment / da;x(e" -I- e'^)(a;)cont summing over respectively occupied, 
Eq. (31), and non-occupied, Eq. (32), states, and take the difference. 
This yields the A{A) defined in Eq. (36), which is expected to be zero in 
the limit A ^ 0. The numerical results for ^(A) for 3 GeV < A < 6 GeV 
are shown in Fig. 4. We see that A{A) tends to zero with increasing A 
but very slowly. (We checked that one obtains within few percent the 
same result from integrating the curves for a;(e" -I- e'*)(x)cont in Fig- 3.) 

This slow vanishing of the "anomaly" (in the sense of [49]) ^(A) is, 
in fact, precisely what we expect. In order to see that, we evaluate 
the theoretical expression (35) for ^(A) in the same way as done in the 
numerics, namely for the Wood-Saxon regulator 
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FIG. 3: The continuum contribution to 
a;(e" + e'*)(a;) as function of x from sums over 
occupied ("occ") and non-occupied ("non") 
levels for different intermediate cutoffs A. 
The arrows indicate the tendency of the 
curves with increasing cutoff. 

A(A) 
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" , 1 , , 








i?(w,A) 



1 



1 + exp[ 



A)/£] 



(39) 



-0.5 



3 4 5 A (GeV) 

FIG. 4: The difference A{K) defined in (33) 
as function of the intermediate cutoff A. It 
must vanish for A — > oo according to (36). 
The thick curves are the numerical results. 
Solid: for "usual" Pauli-Villars masses used 
throughout. Dashed: for a set of smaller 
masses Mi = 389 MeV and A/2 = 700 MeV. 
The thin curves are the respective analytical 
results from Eq. (35) . 



with the box-size dependent width e — 4/i?box (here i?box = 12 fm), and 
integrate in spherical coordinates over k up to |k| < Anum = 8.9 GeV. 
The result is shown in Fig. 4 as dashed line, and we observe again an 
impressive agreement of the analytical and numerical results. 

To test further the theoretical result (35) we repeated the calculation with M\ — 389 MeV and M2 — 700 MeV. With 
smaller (compared to the "usual" ones fixed in the sequence of Eq. (23) and used throughout) Pauli-Villars masses 
one subtracts "more". (The opposite is evident, for M\^2 ^ 00 one recovers unregularized, divergent expressions.) 
Therefore, continuum contributions and consequently ^(A) are smaller. Also in this case the analytical and numerical 
results agree remarkably, see Fig. 4. 

It happens that ^(A) is much larger then the continuum contribution whose computation it hampers, namely (in 
the numerator of the undesired, anomalous terms cancel out) 



/ da;a;(e" -|-e'')(a;,A)c 



/ da;a;(e" e'^)(a;, A)co„t, 



/ da;a;(e" -t-e'^)(2:, A)co„t,occ + / da; x(e« -t- e'i)(a;, A) 



cont . non 



> 



10 for usual Pauli-Villars masses 
. 3 for small Pauli-Villars masses 



(40) 



for 3 GeV < A < 6 GeV. The result (40) indicates that, the better the condition Mi -C A ^ Anum is realized, the 
less the undesired anomalous difference is disturbing. Thus, were we able to use intermediate cutoffs satisfying more 
convincingly the condition A 3> M^, it would be possible to perform reliably the extrapolation A — > 00. However, 
here we are restricted, since a numerical cutoff Anum of substantially more than 10 GeV would require unacceptably 
long computing times. Unfortunately, this means that with justifiable effort one cannot establish the equivalence of 
summations over occupied and non-occupied levels. 

The calculation in Eqs. (32-40) of [22] reveals how this contribution, which only slowly vanishes with increasing A, 
is distributed in x. The anomalous terms appear at a; < (a; > 0) when computing a;(e" -Ve'^){x, A, Mi) from occupied 
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(non-occupied) states. Due to the smearing procedure, however, the effects of the anomalous terms penetrate also in 
the respectively opposite x-regions, and yield the picture in Fig. 3. 

^From this observation it is clear how one can evaluate a;(e" + e'^){x). We have to switch off the smearing, and for 
a; > (a; < 0) we have to sum over occupied (non-occupied) states. 

The observation that for x > (a; < 0) the sums over occupied (non-occupied) states converge faster than in the 
opposite a;-regions was already noticed in [45]. Still, in all examples encountered so far, the convergence in those 
slower-convergence- x-regions was fast enough in order to establish reliably the equivalence of results from summations 
over occupied and non-occupied states [45-47, 49, 50, 53]. Here we face for the first time the situation where this is 
not possible, which is not surprizing as it is also the first time one has to deal with Pauli-Villars-masses of 0(1 GeV), 
and we are forced to give up this strong test of the numerical results. 

The price to pay for giving up the smearing step (26) is less serious. One may turn discontinuous results into final, 
continuous ones, see below. 

Fig. 5a shows the results for the continuum contribution to the antiquark distribution, a;(e" -I- e'*)(a:)cont, non (i-e. 
a;(e" -I- e'')(x)cont at negative x) obtained from the sum over non-occupied states without smearing, for 7 ^ in (26). 
The discountinuous nature of the curve is apparent (we connected the points to guide the eye). It is remarkable that 
our results practically do not depend on the intermediate cutoff A. In other words, we observe a very fast convergence 
in A. It is clear that whatever procedure we use to "smooth" the curve, it will introduce larger numerical uncertainties 
then the extrapolation in A, and therefore refrained from performing the latter. The solid curve in Fig. 5a shows the 
final, smoothened curve obtained from fitting the discontinuous curves. 

Fig. 5b shows the corresponding results for a;(e" -I- e'^)(x)cont, occ obtained from the sum over occupied states. Also 
here we observe a stable convergence in A, except for very small x ~ 0.05, where anyway the finite box method is not 
reliable, and the large- A^c approach not applicable [44, 45]. The solid curve in Fig. 5b shows the final result obtained 
from fitting the discontinuous curves. (Since the continuum contribution to the quark distribution is far smaller than 
to the antiquark distribution, we computed it with a more coarse-grained resolution in x.) 

What about the respective slow-convergence-x-regions? Here, without smearing, one observes fluctuations similar 
to that in Fig. 5, but two orders of magnitude larger, which are due to the anomalous terms, see the discussion above. 

Having switched of smearing, we could have equally computed (e" -I- e'^){x) (since now computing x{e'^ + e'^){x) or 
computing (e" -I- e'^){x) and multiplying by x, of course, commutes). Still it is preferable to calculate x(e" + e'^){x), 
and the reason is evident from Fig. 5. With increasing x the 0-functions in (25) allow to include more discrete states, 
and the fluctuations due to the discrete basis diminish. And vice versa, at small x only few long-wave-length states 
contribute, which enhances the fluctuations (and is the reason why here one is particularly sensitive to details of the 
finite box method, as mentioned above). Weighting the function by x suppresses the fluctuations in the small- a; region 
yielding a less fluctuating curve which can be "smoothened" more reliably, as is demonstrated by Fig. 5. 

The final result for the regular part of (e" -I- e'^){x) is shown in Fig. 5c. For completeness we show how it is 
composed of the contributions from the discrete level. Fig. 1, and the continuum contribution. Figs. 5a and b. (In the 
smoothening step yielding the final (solid) curves in Figs. 5a and b we have build in the constrain that the regular 



part to the continuum contribution behaves as (e" 



const for X ^ 0. 



As can be seen from Fig. 5c, the regular continuum contribution to (e" -I- e'^){x) is small compared to the discrete 
level contribution which dominates the final result. Although we could not check the equivalence of the regular 
continuum results from summations over occupied and non-occupied states, we still were able to clearly demonstrate 
that the numerical calculation is well under control, which gives confidence into the final results. 
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FIG. 5: The regular continuum contributions to (a) x{e'" + e'*)(j:), and (b) x(e" -I- e''){x) as functions of x for the intermediate 
cutoffs A = 3, 4, 5 and 6 GeV, and the respective final, smoothened results, (c) The total result for the regular part of 
(e" -|- e'*)(a;) as function of x, and how it is decomposed of the discrete level and continuum contributions. 
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FIG. 6: (a) The total result for (e" — e'*)(a:) as function of x, and how it is decomposed from the discrete level and continuum 
contributions, (b) A detail on how the final result for the continuum contribution in Fig. 6a comes about: (e" — e'*)(2;, A)cont 
as function of x for different intermediate cutoffs from sums over respectively the occupied and non-occupied states. 

VII. CALCULATION OF (e" - e'^)(x) 

The flavour non-singlet combination (e" — e.'^){x) appears only at subleading order in the l/iVc expansion, when one 
includes "rotational" corrections. This flavour combination is UV-finite and does not need to be regularized. Fig. 6a 
shows the final result for (e" — e'^){x). It is also shown how the total result is decomposed of respectively the discrete 
level and continuum contributions. 

We observe a fast and stable convergence of the continuum contribution (e" — e'^)(a;, A)cont with increasing interme- 
diate cutoff A. This is demonstrated in Fig. 6b where only the region of strongest A-dependence around \x\ < 0.3 is 
shown. Notice that the different curves in Fig. 6b would be nearly indistinguishable on the scale of Fig. 6a. After the 
extrapolation A — > oo one obtains from the sums over respectively over occupied and non-occupied states, Eqs. (17) 
and (18), results which coincide to within an accuracy of about 1%. 

The final result for (e" — e'^){x) shown in Fig. 6 reveals that it is a regular function of x. In particular, no S{x)-type 
singularity appears in this fiavour combination. 

In order to separate consistently different flavours, e"(a;) and e'^{x), it would be necessary to consider also rotational 
corrections to the leading large- A'^c structure (e" + e'^){x). Applying straight-forwardly the procedure described in 
Sec. Ill which lead to the expressions (15-18) one obtains rotational corrections to (e" + e'^){x) consisting only of 
incomplete double sums. This is a general feature encountered whenever one considers 1/iVc-corrections in the model 
to those parton distribution functions which appear already at leading order of the large- iVc expansion [50]. Below, 
when discussing the final results for e°'{x), we shall follow the suggestion [24, 48] to discard such terms. 

VIII. DISCUSSION OF RESULTS FOR e°-{x) AND SUM RULES 

It is interesting to observe that e"'{x) clearly respects the large- A^c predictions for the fiavour dependence, Eq. (7). 
The "large" flavour combination {e^ + e'^) (x) is indeed much larger than the "small" flavour combination (e" — e'^) (x) , 
see Fig. 7a. As a consequence one finds e"(a;) « e^(x) within an accuracy of about 30 %, see Fig. 7b, which is precisely 
what one generically expects from next-to- leading order corrections in an l/A'^j.-expansion with Nc = 3. The large- A'c 
predictions are even more convincingly realized in the case of antiquaries, see Fig. 7c. 

In order to gain some more intuition on the model results for e'^{x) it is instructive to compare them to ff{x) 
computed in the same model [46, 47]. The xQSM results for /f (x) agree to within an accuracy of about 30% with 
parameterizations performed at comparably low scales [55, 56]. 

Figs. 8a-d show xe°'{x) in comparison to xfi(x) for a = u, d, u, d at the low scale of the model. (In this figure it can 
be seen best that in the model parton distribution functions have a non-zero support also for x > \ where, however, 
they are proportional to exp(— const A'ca;) [44]. Since our results refer to the large- A'^c limit, there is conceptually no 
problem. In practice, even for = 3 the distribution functions are negligibly small for a; > 1, see Fig. 8.) 

It is remarkable that the e°'{x) are sizeable only for x < 0.5, in contrast to the fi{x) which extend also to larger 
values of x. In the region olx 0.5 the quark distributions e'^(x) are about half the magnitude of the ff{x). However, 
the antiquark distributions e^{x) and flix) are of comparable magnitude. 

This comparison is interesting also because it was shown that in the non-relativistic limit e'^{x) and ff{x) coincide 
[4, 22]. In the xQSM, which is a relativistic model, one is far away from a non-relativistic scenario, see Figs. 8a and b. 
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FIG. 7: (a) Comparison of (e" +e'*)(a::) and (e" -~e''){x) as functions of x. (No attempt is made to indicate the (5(2;)-contribution 
in the flavour singlet case.) (b) The quark distributions e'"{x) and e''{x). (c) The antiquark distributions e'"{x) and e''{x). 
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FIG. 8: The different flavours of xe°'(x) as functions of x in comparison to from the same model, Refs. [46, 47]. 



Next let us discuss the sum rules (5, 6) virhich were analytically proven to be satisfied in the model in Refs. [22-24]. 
For the first moments we obtain here 



dx (e" + e'^){x) = 1.76iev + 0.11co„t,rcg + (4.88 + 0.196)cont, 



j dx (e" - e'^){x) = 0.303iov + O.OOlcont = 0.30 , 



sing 



= 6.95. 



(41) 
(42) 



where we indicate how the numbers are composed of respectively the discrete level, and the regular and singular 
continuum contributions (the last corrected for the finite box size effects we know, see Sec. V). The results are 
in agreement with phenomenology, see (5, 6). Notice, that the result in (41) is strongly sensitive to details of 
regularization. With the exact result for the coefficient C we obtain J da; (e" + e'*)(a;) ~ 12, see Sec. V. 

While the sum rules (5, 6) are formally satisfied in the model [22-24] and numerically in satisfactory agreement 
with phenomenology one must admit a shortcoming. In QCD the sum rules (5, 6) are saturated solely by the 6{x)- 
contribution, see Sec. II. In contrast to that in the model about 20% of the result in (41) are due to the regular part 
of (e" + e'^){x), while the total result in (42) is due to the regular (and only) part of (e" — e'^)(x). 

Another shortcoming is that in the chiral limit to which our results refer the second moments (4) for (e" ± e'^){x) 
vanish in QCD. Instead, we obtain in the model 



J dx a;(e" + e'^){x) = 0.258^ - 0.024cont = 0.23 , 



da; a;(e" - e'^){x) = 0.089iev - 0.019cont = 0.07 . 



(43) 
(44) 



These two shortcomings have the same origin. 
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Both, the fact that the sum rule (3) is saturated solely by the (5(.T)-contribution and the sum rule (4), follow from 
applying explicitly the QCD equations of motion. At this point the xQSM and actually any model is overburdened. 
Effective model approaches do not respect the QCD equations of motion, rather they satisfy the respective model 
equations of motion. Therefore it is not surprizing to observe such sum rules not to be satisfied literally. Still, one 
may explore model equations of motions and reinterpret the sum rules (3, 4) in the model terminology [22, 24]. 

IX. COMPARISON TO PREVIOUS CALCULATIONS IN THE xQSM 

The first (approximate) calculation of e"'{x) in the xQSM was reported in [22]. There the regular part of (e" + e'^)(a;) 
was approximated by the discrete level contribution, and rotational corrections were discarded. The accuracy of these 
approximations was estimated to be within 30 %. In fact, the (regular) continuum contributions are small compared to 
the respective discrete level contributions, see Figs. 5c and 6a, and the rotational correction (e" — e'')(a;) <S (e"+e'')(x), 
see Sec. VIII. Thus, the estimates of [22] did indeed approximate e"(a;) within the claimed accuracy. Of course, this 
became clear only after exact calculations were presented. 

The first exact calculation of e"(x) in the xQSM was presented in [24], and the present work confirms the results 
obtained there. The main quantitative difference between our work and [24] is that there it was possible to handle 
much larger Pauli-Villars masses (which requires, see Sec. IV, more computing power available in [24]). 

The choice of Pauli-Villars masses is relevant for the continuum contribution {e^ + e'^){x). As a consequence of the 
larger Pauli-Villars masses in [24] an about two times larger coefficient of the (5(x)-contribution was obtained. But 
the results for the regular contribution to e"(a;) obtained here and in [24] practically agree. In fact, these results are 
strongly dominated by the discrete level contributions, see Figs. 5c and 6a, and at this point the slightly different 
value of the constituent mass M ~ 375 MeV used in [24] (compared to M = 350 MeV used here) is more decisive than 
the different choice of Pauli-Villars masses. 

The present work extends in several aspects Ref. [24] , which was up to now the most detailed and complete study 
of e"'{x) in this model. 

• We have shown that the coefficient C of the S{x)-tevm in (e" -|- e'^)(a;) can be determined exactly by means of 

an analytical calculation [22] . 

• We have determined finite box size corrections for the coefficient C of the S{x)-teim in (e" -|- e'^){x). 

• We demonstrated the equivalence of results obtained from summations over occupied and non-occupied states, 
wherever possible. (In [24] the equivalence was demonstrated only for the coefficient C.) 

• Where this was not possible, namely for the regular continuum part of (e" -|- e'^){x), we were able to explain 
why, and to demonstrate that nevertheless the involved numerics is under analytical control. 

In view of the complexity of the task — to deal numerically with a (5(a;)-term, to use a double Pauli-Villars 
subtraction with large Pauli-Villars masses — the present work provides an important supplement to Refs. [22-24] . 

X. COMPARISON TO RESULTS FROM OTHER MODELS 

The distribution function e°'{x) was studied also in other non-perturbative model-approaches, such as bag [1, 18] or 
spcictator [19] models, as well as in 1+1 dimensional toy-models or perturbative one-loop model calculations [20, 21]. 
In Fig. 9 we compare our results to the MIT bag model calculation [1] and the spectator model [19]. 

In the bag model version used in [1] the nucleon is assumed to consist of 3 non-interacting, massless quarks confined 
to the interior of a 3D spherical cavity ( "bag" ) by imposing appropriate covariant boundary conditions which model 
confinement and mimic gluonic effects. The model is relativistic. The flavour dependence is due to the assumed 
SU(2)[iavourXSU(2)gpin spin-flavour-symmctry of the quark wave fimctions such that e"{x) = 2c{x) and e/{x) = c{x), 
and analogously for antiquarks, with e(a;) as introduced in [1]. There is no (5(a;)-contribution in the bag model, and 
since the quarks are massless, the e°'{x) computed in [1] correspond to e^^^{x) in Eq. (2). Interestingly, the pure 
twist-3 ("interaction-dependent") nature of this contribution is refiected in the bag model by the fact that e'^{x) is 
due to bag-boundary effects. The sum rules (3, 4) are not satisfied in the bag model because the QCD equations of 
motion are modified in the bag [1], c.f. the discussion in the xQSM in Sec. VIII. 

The bag gives also rise to antiquark distributions, however, to unphysical ones since the unpolarized antiquark 
distributions in the bag model violate positivity, i.e. in this model ff {x) < is found. Only valence quark distributions 
are considered to be physical [60] . Keeping this mind, we plot also e^{x) from the bag model in Fig. 9c and d. 
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FIG. 9: Comparison of predictions for e"(x) from different models. Solid line: xQSM, computed in this work. Dashed line: 
bag model [1]. Dotted line: spectator model [19]. 



In spectator models parton distribution functions are modelled by introducing a unity in the shape of \n) (n| = 1 
into the definition, here Eq. (1), where {\n)} is a complete set of intermediate states, and modelling this complete set 
of states by e.g. a diquark state. Hereby the diquark state can be taken to be, e.g., in a spin or spin 1 state, which 
is referred to as respectively scalar and vector diquark. Both were considered in [19] which yields, in spite of the 
SU(2)flavourXSU(2)spin spin-flavour-symmetry assumed also there, to a non-trivial flavour dependence in Figs. 9a, b. 
Antiquark distributions can be included by introducing more complex intermediate spectator states, but were not 
considered in [19]. Also the spectator model does not respect the sum rules (3, 4). 

When comparing the different models of e°(x) in Fig. 9 one has to keep in mind that the low scales in the various 
models differ somehow. Bag model results refer to a scale of about 0.4 GeV, spectator model results to about 0.5 GeV, 
and the scale of xQSM results is roughly 0.6 GeV, see Eq. (10). 

The three models agree on that e"(a;) is positive and sizeable for x < 0.5, see Fig. 9a, while e'^{x) appears to be 
smaller (and exhibits in the spectator model even a remarkable zero around x 0.3), see Fig. 9b. The e^{x) are much 
larger in the xQSM model than in the bag model, see Figs. 9c and d. Only in the xQSM there is a ^(a;)-contribution 
at X = (no attempt was made to indicate this singularity in Fig. 9). 



XI. CONCLUSIONS 



A study of the distribution functions e°'{x) in the xQSM was presented which completes previous works [22-24]. 
Two particular features not encountered before in calculations of other parton distribution functions complicate the 
computation of e°'{x) in the xQSM: the appearances of a (5(a;)-singularity and of quadratic UV-divergences whose 
regularization requires a double Pauli-Villars subtraction with a large second Pauli-Villars mass M2 — 0{1 GeV). We 
have demonstrated in detail that in spite of these complications a reliable, controlled numerical calculation is possible. 
Our results confirm qualitatively and quantitatively previous calculations of e°'{x) in the xQSM [22-24]. 

Moreover, our study reveals several interesting results. We have demonstrated that the coefficient C of the S{x)- 
singularity in (e" + e'') (x) can be calculated analytically in the model [22] . As far we are aware this the only case where 
it is possible to compute exactly a (though admittedly unusual) contribution to a parton distribution function in the 
xQSM. This means that the coefficient C of the (5(a;)-term is exactly proportional to the quark vacuum condensate, 
and thus ultimately connected to chiral symmetry breaking [22] . 

Another interesting, but more technical byproduct of our study is that we have been able to quantify exactly the 
finite-box-size corrections to a quantity in the model, namely to the coefficient C (which are of the order of few %). 
This is also, to best of our knowledge, unique in the model. 

The xQSM predicts e^{x) to be positive and sizeable — reaching half the size of /"(x) — in the region of a; < 0.5, 
and e'^{x) somehow smaller. Remarkably, in this region of x the e^{x) appear to be of similar magnitude as the fi{x). 
Our predictions for the quark distributions e'^{x) are in rough qualitative agreement with results from bag [1, 18] or 
spectator [19] models in which, however, no 5(a:)-contribution appears. 

Here, following the practical point of view of Ref. [1], we computed e°'{x) benefiting from a special property of this 
twist-3 quantity, which allows to define it in terms of quark fields only, and thus makes possible studies in models 
without explicit gluon degrees of freedom [1]. However, by means of QCD equations of motion e"'{x) can be reexpressed 
such that they explicitly depend on gluon fields. The question whether quark models, like the xQSM, are nevertheless 
able to provide useful estimates for e°'{x), at least in certain regions of x, can be clarified only by experiment. 
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